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Abstract 

A new method to measure the angle of the CKM unitarity triangle using a Dalitz analysis of 
the three-body decay of the D° meson from the B — > D°K process has been introduced recently 
The method employs the interference between D° and D° to extract both the weak and strong 
phases. We present a study of statistical and systematic errors of this method using Monte Carlo 
simulation. A fit of the D° Dalitz distribution has been performed based on a 140 fb _1 sample 
collected by the Belle detector. We observe the indication of a CP asymmetry with the significance 
of 2.4 standard deviations. The 90% confidence interval obtained for the ^3 value is 61° < ^3 < 
142°. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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INTRODUCTION 



The determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements [l| is im- 
portant to check the consistency of the Standard Model and search for new physics. Various 
methods using B — > DK decays have been introduced to measure the unitarity trian- 
gle angle 03 but these require knowledge of the branching fractions for different B — > DK 
modes and therefore suffer from systematic errors due to detector response and background 
contamination. 

A novel technique based on the analysis of the three-body decay of the D° meson [3] 
does not involve absolute branching fraction measurements and is, therefore, free from cor- 
responding detector systematics. At the same time, the uncertainties in this method due to 
the D° decay model can be controlled. 

This method is based on two key observations: neutral D° and D° mesons can decay to 
a common final state such as Ks^ + tt~, and the decay B — > D°K^*> can produce neutral 
D° mesons of both flavors via b — > cus and b — > ucs transitions (see Figure 1), where the 
relative phase between the two interfering amplitudes is the sum 5 + 03 of strong and weak 
interaction phases. In the charge conjugate mode, the relative phase is the difference 5 — 03 
of these phases, so both phases can be extracted from the measurements of such B decays 
and their charge conjugate modes. The annihilation diagram contribution is neglected here 
since it affects only the extracted value of the strong phase 5. In our case, the measurement 
is based on the analysis of Dalitz distribution of the three body final state of the D^*> meson. 

In this paper, we will discuss in detail the decay B + — > D°K + state only, although a 
similar approach can be applied to other processes such as B° — > D°K*° or B^ — > D*°K ± . 



DESCRIPTION OF THE METHOD 

In the Wolfenstein parameterization of the CKM matrix, the amplitudes of the two 
diagrams shown in Fig. 1 that contribute to the decay B + — ► D°(D°)K + are given by 
Mi ~_ VcbV*, ~ AX 3 and M 2 ~ V ub V* s ~ AX 3 (p + irf). These amplitudes interfere if the D° 
and D° mesons decay into the identical final state K s 7t + tt~; we denote the admixed state 
as D°. Assuming no CP asymmetry in D° decays, the amplitude of the B + decay can be 
written as 

M+ = f(m 2 + , m 2 _) + ae l ^ +l5 f{m 2 _ , m\), (1) 

where m 2 + and m 2 _ are the squared invariant masses of the K s 7r + and K s ir~ combinations, 
respectively, and / is the complex amplitude of the decay D° — > K s -k + -k~ . The absolute 
value a of the ratio between the two interfering amplitudes is given by 

a = j^ u& ffi • M ~ 0.08 ■ 1.04/0.2196 ■ 0.35, 

\VcbV* s \ [ai| 

where the color suppression factor IfWail can be estimated from the Belle measurements of 
the color suppressed B° — > D°K° and color allowed B~ — > D°K~ Q decays: 



\a 2 \ / Br(B° -> D°K°) 

H ~ V Br ( B ~ D ° K ~) ~ °' 35 ' 

Though there are estimations of the a value as large as 0.2 |7(, in our Monte Carlo (MC) 
simulations, we have used the conservative value of a = 0.125. 
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Similarly, the amplitude of the charge conjugate B decay is 

M_ = f{m 2 _,m 2 + ) + ae- i ^ +i5 f{m 2 + ,m 2 _). 

Once the functional form of / is fixed by a choice of a model for D° — > K s tt + ti~ decays, 
the Dalitz distributions for B + and B~ decays can be fitted simultaneously by the above 
expressions for M + and M_, with a, 03, and 5 as free parameters. 

Ref. ||| suggests a model-independent way of determining 3 via a binned Dalitz plot 
analysis. However, this method is likely to introduce higher statistical errors in our present 
statistics-limited stage, and certain technical issues involving backgrounds and reconstruc- 
tion efficiency have to be solved. Here, we use a model-dependent approach with an unbinned 
maximum likelihood fit of the Dalitz distribution, making optimal use of our small event 
sample. 

STUDY OF SENSITIVITY AND MODEL UNCERTAINTIES 

The accuracy of the 03 determination and possible systematic effects are determined from 
a MC simulation, where we generate a large number of D° — ► K s tt + tt~ decays according 
to a selected model of the D° decay, then perform the unbinned likelihood fit of the Dalitz 
distribution to obtain the parameter values and their errors. 

The Dalitz plot for the decay B + — ► [K s -k + 'k^]^ K + in the variables m 2 + and m 2 _ was 
generated according to relation (JTJ), with the D° — > K s 7i + ti~ decay amplitude / described 
by the set of resonances listed in Table |1J The model used is similar to that in the CLEO 
measurement [H except that tensors and If* (1680) are excluded. The amplitude / in this 
model is represented by a coherent sum of N two-body decay amplitudes plus one non- 
resonant decay amplitude: 

N 

f{m 2 + ,m 2 _) = J2 a j eiaj M m l> m -)+ beil3 , 

3=1 

where N is the number of resonances, Aj(m 2 i _,m 2 _), aj and otj are the matrix element, 
amplitude and phase, respectively, for the j-th resonance, and b and (3 are the amplitude 
and phase for the non- resonant component. The total phase and amplitude are arbitrary. 
To be consistent with CLEO, we have chosen the K s p mode to have unit amplitude and 
zero relative phase. The description of the matrix elements A, follows Ref. 

The simulated Dalitz _plot of the D° decay is shown in Fig. El The corresponding distri- 
bution for the admixed D° decay contains the information about the interference phase and 
amplitude between D° and D° that can be extracted from the fit to the Dalitz distribution. 

Though in general the whole Dalitz plot is sensitive to the variations of the phase between 
D° and D°, the contribution from the interference is particularly high in certain regions. 
Since the maximum of the D° plot density is in the region of K*(892) + , the most sensitive 
area is the region of double Cabibbo suppressed if* (892) _ , which interferes with the maxi- 
mum amplitude of color suppressed D°. Another sensitive part is the region of K s p. Fig. 0] 
demonstrates the CP asymmetry on the most sensitive regions of the Dalitz distribution. 

For the estimation of the statistical error of 3 , samples of 10 4 events were generated 
with total phase 9 varying from to 27r (8 = 5 + 03 for B + — > D°K + and 9 = 5 — 03 for 
B~ — > D°K~ decay) and the resulting D° distributions were fitted with a and 6 as free 
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parameters. This is equivalent to the case where both B + and B decays are fitted with 
(a+, 9+) and (a_, 0_) as free parameters and the weak phase 03 is extracted as 




The error of 03 is slightly overestimated in this case since we do not constrain the relative 
amplitudes a + and a_ to be equal. 
The error of 03 was determined as 

aj 3 (0 3 , 6) = -y e {8 + 3 ) + oi{6 - 3 )), (2) 

where <jg is the error of the parameter 9 extracted from the fit. The corresponding statistical 
error <t^ 3 obtained with relation (J2J) ranges from 3.0° to 3.7° with a mean value of 3.4° 
(averaged over all 3 and 5). 

We expect about 10 3 B + — ► [K s -k + 'k^] i ^oK + decays from a 1000 fb _1 sample will result in 
a 10° statistical accuracy in 3 . The statistics of B decays collected at the present time by 
5-factories is therefore insufficient to make a 03 measurement with a reasonable precision. 
However, the 140 fb _1 sample of Belle could provide the evidence of a CP asymmetry. 

We expect the model used for the D° — > K s n + n^ decay to be one of the main sources 
of systematic errors of 03. We estimate this contribution by fitting our simulated samples, 
described earlier, to five alternate models for the decay amplitude f(m+, m 2 _) (see Table ITT|). 
It can be seen that while Blatt-Weisskopf form factors and q 2 resonance width dependence 
have a rather small effect on the model error, the exclusion of wide resonances and doubly 
Cabibbo suppressed if* (892) causes a large bias of 3 . We take the value of 10° as a rough 
estimate of the model uncertainty for the method. This value was obtained by taking only 
the narrow resonances (including doubly Cabibbo suppressed if* (892)) and approximating 
the others with a flat non-resonant term. This is motivated by the fact that the amplitude 
of the narrow resonances is well approximated with a Breit-Wigner shape, whereas wide 
states (r ~ M) contain substantial theoretical uncertainties. 

Though doubly Cabibbo suppressed if* (892) has a small amplitude compared to other 
resonances (its fit fraction is 0.34% jij), the error in its amplitude transforms into a large 
systematic error since the whole area of the resonance interferes with the Cabibbo allowed 
if* (892) component of the color suppressed amplitude. We have investigated this influence 
by generating a sample of 10 4 events according to the model of Table |l] with subsequent 
fitting of this Dalitz plot with the same model where the amplitude of doubly Cabibbo 
suppressed if*(892) is changed. This analysis indicates that the statistical error of 0.02 on 
the if*(892) amplitude obtained by CLEO transforms into approximately 3° bias in 03. 
The errors of other amplitudes are expected to contribute even less, since they lie mainly 
outside the maximum sensitivity region of the Dalitz distribution. 

Another source of systematics comes from various kinds of background. We expect three 
main background contributions to affect the determination of 03: 

• Background uniformly distributed over the Dalitz plot (combinatorial background), 

• Background from B^ — > D Q {D°)'K ± events with a pion misidentified as a kaon (pion 
background) , 

• Background from 5 ± — ► D°(D°)K ± events with a wrongly associated kaon of opposite 
sign ("wrong tag" background). 
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To consider the worst case scenario, our MC samples were generated with one or another 
of these backgrounds and then fit without any background component. Neglected combina- 
torial background at the level of 5% results in an 8° bias in 3 . Neglected Dn background, 
even at the level of 20%, does not induce an appreciable bias in 3 . The contribution of 
wrongly associated kaons, which tag the flavor of D°, has a larger effect on systematics, since 
half of such kaons have the wrong sign and such events will be misinterpreted as decays of a 
D meson of opposite flavor. These D 0, s will distort the most sensitive region of the Dalitz 
plot. According to our simulation, the 5% admixture of wrong tag background introduces a 
12° bias of <p 3 . 



DETERMINATION OF D° -> K s ttit DECAY MODEL 

For our MC simulation we have used the D° decay model measured by CLEO based on 
their data sample of 5300 events However, the integrated luminosity collected by Belle 
can offer a much bigger sample of D° events. We have analyzed 78 fb _1 of BELLE data 
looking for decays of D*(2010) + — > D°tt + , D° — > K s 7t + tt~ (and its charge conjugate). The 
selection was based on the D° candidate invariant mass (\M Ks7T7r - M D o\ < 9 MeV/c 2 ), and 
the mass difference of the D* + and D° candidates (4.39 < AM < 6.19 MeV/c 2 ). We require 
the CM momentum of the Dtx system to be above 2.7 GeV/c to avoid contamination from 
combinatorics of bb decays. We have selected 57800 events satisfying these requirements. 
The fit of the distribution yields a background fraction of 5.6%. 

The Dalitz plot for the selected D* — > Dtt — > (KsTi'Ti')^ candidates and its projections 
are shown in Fig. EJ The data were fitted with a list of resonances from the CLEO model 
plus two scalars (<Ji and 02) in the tttt channel with free mass and width. We assumed 
the background and efficiency to be uniform over the allowed phase space. Our fit result is 
shown in Table ITTT1 



EVENT SELECTION 

The Belle detector has been described in detail elsewhere It is a large-solid-angle 

magnetic spectrometer consisting of a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC) for charged particle tracking and specific ionization measure- 
ment (dE/dx), an array of aerogel threshold Cerenkov counters (ACC), time-of-flight scintil- 
lation counters (TOF), and an array of 8736 CsI(Tl) crystals for electromagnetic calorimetry 
(ECL) located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron flux return located outside the coil is instrumented to detect Kl mesons and identify 
muons (KLM). 

Separation of kaons and pions is accomplished by combining the responses of the ACC 
and the TOF with dE/dx measurement in the CDC to form a likelihood C(h) where h is a 
pion or a kaon. Charged particles are identified as pions or kaons using the likelihood ratio 
(PID): 

PID < A '» - zmtw) 1 PID w = zmtW) = 1 " ™- 

The data sample of 152 • 10 6 BB events (140 fb _1 ) collected by the Belle detector was 
processed to search for B ± — > D°K ± decays with D° decaying into the K s tt + tt~ final state. 
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In addition, we select — > D ^ and — > D*(2007)°7r ± decays which serve as test 
samples for our fitting procedure and are used for background analysis. 

Neutral kaons are reconstructed from pairs of oppositely charged pions. We require the 
reconstructed vertex distance from the interaction point in the plane transverse to the beam 
axis to be more than 1 mm and the invariant mass M nn to be \M 7Tn — Mk s \ < 10 MeV/c 2 . 

Selection of B candidates is based on the CM energy difference AE = E beam and 

B meson beam-constrained mass M bc = ^ E 2 eam — (^Pi) 2 , where E beam = m T /2 is the CM 
beam energy, and E^ and pi are the CM energies and momenta of the B candidate decay 
products. We select events with M bc > 5.2 GeV/c 2 and \AE\ < 0.2 GeV for the analysis. 
The cuts for signal candidates are 5.272 < M bc < 5.288 GeV/c 2 and \AE\ < 0.022 GeV. In 
addition, we apply a cut on the invariant mass of the D° candidate: |M^ sVr7r — m£>| < 11 
MeV/c 2 . 

To suppress the continuum background, we require | cos 6*4^,. | < 0.8, where t hr is the 
angle between the thrust axis of the B candidate and the rest of the event. For additional 



background rejection, we use the Fisher discriminant composed of 11 parameters [ll|: the 
production angle of the B candidate, the angle of the B thrust axis relative to the beam 
axis and nine parameters representing the momentum flow in the event relative to the B 
thrust axis in the CM frame. 

The selection efficiency for these cuts is 12%, based on MC simulation, implying an 
expectation of 158 ± 33 signal events in our data sample. The distributions for AE, M bc 
and Md for the experimental data are shown in Fig. E| The number of B — > D°K events 
satisfying all selection criteria is 138. 

To fit the M bc distribution, we use the empirical background function suggested by the 
ARGUS collaboration 12[ plus a Gaussian signal. The AE distribution is fit with a linear 
background plus two Gaussians representing the signal (with a mean value around zero) and 
background from B — > D ir with a pion misidentified as a kaon (with AE ~50 MeV). For 
the Md fit we use a constant background plus a combination of two Gaussians for the signal. 
The resolutions of the parameters were extracted from the fit of B — > D tt data, described 
below. 

The fit of the AE distribution yields 107 ±12 signal events and 33 ±3 background events, 
consistent with the expectation of 158 ± 33 signal events. The Dalitz plots of the selected 
candidates and their projections onto + , M\ _ and M 2 +vr _ axes are shown separately 
for B~ and B + decays in Fig. Eland Fig. [TO respectively. Note that for B + decays the axes 
are chosen so that the K* band is aligned along the ?/-axis, as for B~ decays. 

The B^ — > D ^ selection cuts are the same as those for B — > D°K, except that the 
cut on K/tc PID probability is modified to select a pion (Pr(K/-n-) > 0.5). The same fit 
functions are used to fit AE, M bc and Md distributions as for B — > D°K (with only one 
Gaussian peak in AE), and the fitted distributions are shown in Fig. |3 The extracted 
resolutions of the selection parameters, used in the B — > D°K fits, are 13 MeV for AE, 3.2 
MeV/c 2 for M bc and 6.3 MeV/c 2 for M D - The Dalitz distribution of the D° -> K s tt + 7T' 
decay from B — > D°7r process is shown in Fig. [7| 



BACKGROUNDS 

To fit the Dalitz plot correctly, backgrounds have to be taken into account in the fit 
model. For each of the background contributions, we need to obtain not only its fraction in 
the event sample, but also the Dalitz plot shape. 



The largest contribution for the decay of our interest comes from the continuum qq 
background (q = u,d,s,c). It includes the background with purely combinatorial tracks, 
and continuum D° mesons combined with a random kaon. This type of background is 
analyzed using the off- resonant statistics collected below T(AS) as well as on-resonance data 
with cos 9thr an d Fisher cuts applied to select continuum. 

The background coming from bb events originates either from B — > D°K decay with 
some final state particles interchanged with combinatorial ones from the decay of the other 
B meson, or from other charged or neutral B decays (possibly with misidentified or lost 
particles). We subdivide the bb background into four types. 

The B — > D°7r process with a pion misidentified as a kaon is suppressed by the require- 
ment on the Kj-ri identification probability and CM energy difference. The fraction of this 
background is obtained by fitting the AE distribution; the Dalitz plot distribution is that 
of D° without the opposite flavor admixture. 

Other decays of charged and neutral B mesons contributing to the background are in- 
vestigated with a generic MC sample. Most surviving candidates are due to decay of one of 
the B mesons to the D*°X state, with some particles taken from the other B decay. 

Events where one of the D° decay products is substituted by a combinatorial kaon or 
pion were studied using the MC data set where one of the B mesons from the T(AS) decays 
into the D°K state. 

Events where a valid D° is combined with a random flavor-tagging kaon are of importance, 
because half of the kaons have the wrong sign: such events will be misinterpreted as decays of 
D°, thus introducing distortion in the most sensitive area of the Dalitz plot. This background 
was studied with the same MC sample as for D° combinatorics. No events of this kind were 
found, which allows us to put an upper limit of 0.4% (at 95% C. L.) on the fraction of this 
contribution. 

The Dalitz plot shapes of all the backgrounds were parameterized with empirical functions 
that included polynomial terms as well as Breit-Wigner peaks. The investigation of the AE 
distribution for all considered backgrounds allows us to conclude that none of them except 
combinatorics in D° decay is peaked in AE. Therefore, a AE fit gives a good estimate of the 
total background rate, which is 23%. The fractions of each contribution are as follows: 14% 
for continuum combinatorics, 3% for D° from continuum, 0.5% for B — > D°ir, 4% for decays 
of B^ other than B — > D°K/n, 2% for decays of neutral B and 0.4% for combinatorics in 
D° decay. 

Another approach to take into account backgrounds is to use events from D° sidebands. 
Since all the backgrounds except D° from continuum and D° decay combinatorics are not 
peaked in the D° invariant mass distribution, the D° sidebands are a good approximation of 
the total background shape. The D° combinatoric background has a small fraction and can 
be neglected, and the shape of the continuum D° background is known and can be added 
separately. However, the statistics in B — > D°K decays is not sufficient to fit the background 
shape with good precision, so we use the B — > D°n sample for this purpose. This approach 
is used together with the flat background distribution to estimate the systematic uncertainty 
caused by our background parameterization. 

FIT OF D° DALITZ PLOT 

We use the unbinned likelihood technique to fit the Dalitz plot of both MC simulated and 
experimental samples. In the fit function of experimental data, finite momentum resolution 
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and efficiency were taken into account. 

The value which is minimized in the fit is the inverse logarithm of the unbinned likelihood 
function 



-21ogL 



log p(m 2 i+ ,m 2 



log 



p(m+, w?J)dm 2 + dm 2 _ 



D 



where m 2 + , m 2 _ are measured (or generated) Dalitz plot variables and p(m 2 +1 m 2 _) is the 
Dalitz plot density being fitted. The Dalitz plot density is a sum of the signal D° — D° 
interference \f(m 2 h ,m 2 _) + ae~ 1 ^ f{m 2 _ ) m 2 + )\ 2 and the background distribution, convoluted 
with a momentum resolution function and multiplied by the efficiency shape. The efficiency 
and resolution functions were obtained with a MC simulation with D° decaying according 
to the phase space distribution. 

To check the consistency of the 03 fit, the same fitting procedure was applied to B — > 
D^°% and D* — > D°n test samples as well as the B — > D°K signal. The idea is to 
check the absence of opposite flavor amplitude a that can appear if the Dalitz plot shape 
is not described well by the fit model. In the case of D* — > D°tt data, we do not expect 
any contribution from opposite flavor (provided flavor tagging is correct), but this sample 
cannot serve as a good fit test since the same data are used to fit the D° decay model. 
So, these data can only indicate the defects of D° decay parameterization, which can affect 
the 03 fit. The sample of D° from B — > D^°ti decay allows for a more reliable test of 
the fit procedure, since it is independent of the sample used to determine the model, and 
the background sources are expected to be similar to those for our signal. In the case of 
B — > D^°n, however, a small opposite flavor amplitude is expected (~ |Ki&Kdl/l^*KTczl) 
with a ~ 0.01 - 0.02. 

For our 03 fit check, the free parameters for Dalitz plot fit were the relative amplitude a 
and total phase 9. 

In the D* — > D°tt data fit, we have used the flat background with the same fraction as 
used in the case of determination of the D° model. The result of the fit is the following: 
a = 0.0041 ± 0.0050, 9 = 123° ± 62° for D° decay and a = 0.0047 ± 0.0050, 9 = 114° ± 52° 
for D° decay, which is consistent with zero. 

The fit results for B -»■ D°n are a_ = 0.067 ± 0.027, 9_ = 228° ± 23° for D° from 
B~ D°n- decay and a+ = 0.065 ±0.029, 9+ = 232° ±24° for D° from B + -> D°n + decay. 
It should be noted that since the value of a is positive definite, the error of this parameter 
does not serve as a good measure of the a = hypothesis. To demonstrate the deviation of 
the opposite flavor amplitude from zero, the real and imaginary parts of the complex relative 
amplitude aexp(— i6) are more suitable. Fig.|H]shows complex relative amplitude constraints 
for B~ and B + data separately. It can be seen from the plot that both amplitudes differ 
from the expected value by more than two standard deviations. 

As another check of the consistency of our fit result with the a = hypothesis, we have 
performed a toy MC study. A Dalitz plot of 1700 events (which corresponds to full statistics 
of B — > D°tt with both flavors included) was generated in 566 trials with the model where 
a = 0, and the fit of these data was performed with a and allowed to float. The distribution 
of the parameter a for these trials is shown in Fig. El It can be seen from the plot that only 
~ 1% of the a values are above 0.06. The bias observed needs further investigation. One of 
the possible reasons can be the variation of the background shape that is not parameterized 
by the fit function. Also, one cannot discount statistical fluctuation. 



10 



Another test sample of D 0, s with the small opposite flavor admixture is provided by the 
B — > D*°7T decay with D*° decaying into D°7i°. The constraints on the complex amplitude 
are shown in Fig. |H] (right). In this case, we do not observe any significant opposite flavor 
amplitude. The result of the fit is a_ = 0.057±0.054, 0_ _= 340±65 for D° from B~ -> D*°n- 
decay and a + = 0.041 ± 0.069, 9+ = 164 ± 100 for D° from B + -> D*°tt + . Though this 
sample is smaller (350 events) and it can not serve a reliable test of fit bias, we observe 
phases different from B — > D°7r case, which can indicate that the fitting procedure itself 
does not introduce a bias to the phase. 

For D° — D° interference fit in B — > D°K, we have performed a simultaneous fit with free 
parameters a, 03 and 5. The result of the fit overlaid on the experimental Dalitz distribution 
projections in shown in Fig. El an d Fig. ^2 for B~ and B + decays, respectively. The fit 
yields a = 0.33 ± 0.10, 3 = 92° t™, 5 = 165° It should be noted that the fit errors are 
quite non-parabolic, since the phase error obtained from the fit is inversely proportional to 
the value of a. Therefore, the errors quoted can not reliably represent the accuracy of the 
measurement. The constraint plots on the pairs of parameters (03, 5) and (a, 3 ) shown in 
Fig. give more information about the accuracy and the statistical significance of the fit. 
The plots show one, two and three standard deviation contours of each pair of parameters. 
The central values are marked with crosses. For each of the plots, the third parameter was 
allowed to float. 

Note the two-fold ambiguity of this method: (03, 5) and (03 + ir, 5 + ir) cannot be distin- 
guished, as they lead to the same total phases. Since there are no reliable measurements or 
theoretical predictions of the value of the strong phase, we keep both values as solutions for 
our data. 

The 3 statistical significance — 2\og L/ L max is presented in Fig. 1131 (left). The a and S 
parameters were allowed to float. 

To illustrate the observation of CP asymmetry and its significance, we have performed a 
separate fit of the B + and B~ data sets. The results of the fits are shown in Fig. El (right) 
as constraints on the complex relative amplitude aexp(i8), as in the case of test sample fits. 
The relative amplitude of the B — > D°n test sample sample fit is also shown for comparison. 
The asymmetry of three standard deviations is apparent. 

The influence of the different sources of systematics on our fit result are presented in 
Table II VI The systematic error is dominated by the D° model uncertainty. The component 
related to the background parameterization was estimated by extracting the background 
shape from the Mo sidebands and by using the flat background distribution. As mentioned 
above, the efficiency shape and resolution were extracted from MC simulation. To esti- 
mate their contribution to the systematic error, we have used the flat efficiency and "ideal" 
resolution, respectively. 

The nonzero opposite flavor amplitude we observe in the B — > D°tt test sample may 
indicate some systematic effect such as background structure or a deficiency of the D° 
decay model. Since the source of this bias is indeterminate, we conservatively treat it as an 
additional systematic effect. 

Since the relative amplitude a is positive definite, its fit value is biased towards the higher 
values by approximately one standard deviation. Consequently, the strong dependence of 
the 03 error on a can cause an overestimation of the 03 determination accuracy. However, 
the constraint plot of a vs. 03 (Fig. ^1 right) shows no significant correlation between the 
parameters. Therefore, we have calculated the 03 and 5 for the value of a = 0.23, which is 
one standard deviation lower than the value obtained in the fit. With model and systematic 
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errors taken into account, and a bias correction applied, our fit result can be presented as: 

• Weak phase 3 = 95° ±™ ± 13° ± 10°, 

• Strong phase 5 = 162° tf 5 ± 12° ± 24°. 

The first error is statistical, the second is experimental systematics and the third is model 
uncertainty. Fig. fT4l shows 90% CL regions for the pairs of parameters (03, 5) and (03, a) 
with systematic uncertainty included. The 90% confidence intervals for the fit parameters 
(including systematic uncertainty) are: 0.15 < a < 0.50, 61° < 03 < 142° and 104° < 5 < 
214°. 

CONCLUSION 

We have applied a new method to measure the CP-sensitive angle 03 using the Dalitz 
analysis of the three-body D° decay in the process B — > D°K. The method is directly 
sensitive to the value of 03 and is free from discrete ambiguities (for the 03 values in the 
range from to 180°). The statistical and systematic uncertainties of this method were 
studied using MC simulation. The statistics of 1000 fb _1 should permit us to determine 03 
with a 10° statistical uncertainty. We estimate the model uncertainty of the method to be of 
the order of 10°. This uncertainty can be significantly reduced with the increase of associated 
experimental data, and practically eliminated by using the data from the cr-factory 0]. 

The first measurement of 03 using this technique was performed based on 140 fb _1 statis- 
tics collected by the Belle detector. The significance of the direct CP violation effect (includ- 
ing systematic effects) is 2.4 standard deviations. The values of 03 and relative amplitude a 
obtained are in agreement with the SM expectation. The 90% confidence interval obtained 
for the 03 value is 61° < 03 < 142°. The accuracy of the measurement can be improved in 
the future by adding other suitable decay modes into the analysis. 
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TABLE I: List of resonances used for D° — > K s tt + tt decay simulation. 



Resonance 


Amplitude 


Phase, (deg) 


K*(892)+^- 


1.56 


150 


K sP ° 


= 1 


= 


K*(892)-ir + 


0.11 


321 


K s uj 


0.037 


114 


K s j b(980) 


0.34 


188 


^ s /o(1370) 


1.8 


85 


^(1430)+7r- 


2.0 


3 


K s 7r + 7r~ non-resonant 


1.1 


160 



TABLE II: Estimation of model systematics. 



Fit model 


(A4> 3 ) max (deg) 


iT* (892), p, nonres 


29.3 


K*(892), p, DCS iT* (892), nonres 


6.7 


K*(892), p, DCS K*(892), /°(980), nonres 


9.9 


= = 1 


3.1 


T(g 2 ) = Const 


4.7 



TABLE III: Fit results for D° — > K s tt + tt decay. Errors are statistical only. 



Resonance 


Amplitude 


Phase, (deg) 


i^*(892)-vr+ 






1.706 ±0.015 


138.0 ±0.9 








1.0 (fixed) 


(fixed) 


i^*(892)+7r- 






(13.6 ±0.8) x 10~ 2 


330 ±3 








(32.8 ± 1.8) x 10~ 3 


114 ±3 


^ s /o(980) 






0.385 ±0.011 


214.2 ±2.3 


K s / (1370) 






0.49 ± 0.04 


311 ±6 


^ s / 2 (1270) 






1.66 ±0.05 


341.3 ±2.3 


^ *(1430)-7r + 






2.09 ± 0.05 


353.6 ± 1.8 


i^ 2 *(1430)-vr+ 






1.20 ±0.05 


316.9 ±2.1 


i^*(1680)"vr + 






1.62 ±0.24 


84 ± 10 


K s a x {M ai = 


535 ± 6 MeV, F ai = 


= 460 ± 15 MeV) 


1.66 ±0.09 


217.3 ± 1.4 


K s a 2 (M a2 = 


1063 ± 7 MeV, T a2 


= 101 ± 12 MeV) 


0.31 ± 0.04 


257 ± 11 


non-resonant 






6.51 ±0.22 


149.0 ± 1.6 
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TABLE IV: Estimation of systematic uncertainties. 



Source 


Aa 


A^3, (deg) 


AS, (deg) 


D° decay model 


0.03 


10 


24 


Background shape 


0.015 


3.2 


2.3 


Efficiency shape 


0.010 


2.4 


0.8 


Momentum resolution 


0.010 


2.5 


0.6 


Test sample bias 


0.02 


12 


12 


Total 


0.04 


16 


27 
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FIG. 1: Diagrams of B+ -» L>°^+ (left) and B+ -► L> if+ (right) decays. 
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FIG. 2: Simulated Dalitz plot of D° — > K s ir + ir decay. 
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FIG. 4: Simulated Dalitz plot of D decay and its slices showing the effect of CP asymmetry in 
B — > D°K decay with the parameters a = 0.125, 03 = 70° and 5 = 0. The two curves represent 
the Dalitz plot density for D° decay from B + — > D°K + (solid line) and B~ — > D°K~ (dashed 
line) processes. 
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FIG. 5: AE, M bc and M D distributions for B -> D°.K" decay. Histo gram represents the data, 
solid line is the fit result, dashed lines show the cut limits for the corresponding distributions. 
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FIG. 6: AE, M bc and Mp distributions for B -> D°vr decay. 
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FIG. 7: Dalitz distribution in variables m 2 K n+ , m? K 
process. 



^_ of D° — > K s 7r + ir decay from B — > D°7r 
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FIG. 8: Constraint plots of complex relative amplitude ae ld for B — > L>°7r (left) and -B — > D*°ir 
(right) decays. Contours indicate integer multiples of the standard deviation. 
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FIG. 9: Distribution of the fitted parameter a in 566 trials of 1700 events each, generated with 
a = 0. 
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FIG. 10: M 2 Kn _, M 2 Kit+ , M 2 +jt _ distributions and Dalitz plot of D° -> K s -k + tt~ decay from 
B~ — > D°K~ process. Thick line represents the fit result, thin line is a = case. 
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FIG. 11: M^ n+ , M^_, M^ +n _ distributions and Dalitz plot of D° -> K s 7v + tt~ decay from 
B + — > D°K + process. Thick line represents the fit result, thin line is a = case. 
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FIG. 12: Constraints on 03 and 5 (left) and a and 03 (right) for all parameters free in the fit. 
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FIG. 14: 90% CL regions for fa and S (left) and a and ^3 (right) with systematic uncertainty 
included. 
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